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We perform a quantum theoretical calculation of the noise power spectrum for a phase measure- 
ment of the light output from a coherently driven optical cavity with a freely moving rear mirror. 
We examine how the noise resulting from the quantum back action appears among the various con- 
tributions from other noise sources. We do not assume an ideal (homodyne) phase measurement, 
but rather consider phase modulation detection, which we show has a different shot noise level. 
We also take into account the effects of thermal damping of the mirror, losses within the cavity, 
and classical laser noise. We relate our theoretical results to experimental parameters, so as to 
make direct comparisons with current experiments simple. We also show that in this situation, the 
standard Brownian motion master equation is inadequate for describing the thermal damping of the 
mirror, as it produces a spurious term in the steady-state phase fluctuation spectrum. The corrected 
Brownian motion master equation [L. Diosi, Europhys. Lett. 22, 1 (1993)] rectifies this inadequacy. 
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I. INTRODUCTION 



Interferometers provide a very sensitive method for de- 
tecting small changes in the position of a mirror. This has 
been analysed extensively in the context of gravitational 
wave detection [0-|| and atomic force microscopes [0,1). 
A key limit to the sensitivity of such position detectors 
comes from the Heisenberg uncertainty principle. The re- 
duction in the uncertainty of the position resulting from 
the measurement is accompanied by an increase in the 
uncertainty in momentum. This uncertainty is then fed 
back into the position by the dynamics of the object be- 
ing measured. This is called the quantum back action of 
the measurement, and the limit to sensitivity so imposed 
is referred to as the standard quantum limit. One of the 
pioneers in this field has been Braginsky in various stud- 
ies of measurement aspects of the fluctuations of light 
caused by the moving mirror ||9|,|l0|. 

In real devices which have been constructed so far, the 
quantum back action noise in the measurement record 
is usually small compared to that arising from classical 



sources of noise. However, as the sensitivity of such de- 
vices increases it is expected that we will eventually ob- 
tain displacement sensors that are quantum limited. The 
quantum back-action noise has not yet been seen exper- 
imentally for macroscopic devices, so seeing it is a topic 
of current interest. Once the standard quantum limit 
has been achieved, this will not be the end of the story, 
however. Various authors have shown that it is possible 
to use contractive states Q, or squeezed light [| 12 1, to 
reduce the quantum back action and therefore increase 
the sensitivity of the measurement even further. 

The interferometer we consider here for measuring po- 
sition consists essentially of a cavity where one of mirrors 
is free to move. This system is also of interest from the 
point of view of cavity QED. Usually cavity QED exper- 
iments require optical cavities where the atomic excita- 
tions and photons in the optical modes become entan- 
gled. The dynamics follows from the interplay between 
these quantum variables. However, a challenging realm 
for cavity QED experiments involves instead an em,pty 
cavity (that is, a cavity containing no atoms or optical 
media) where the photons in the cavity mode interact 
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with the motion of one of the cavity mirrors. In this 
scheme, the position of at least one mirror in the optical 
resonator is a dynamic variable. The coupling between 
the photons and the mirror position is simply the radi- 
ation pressure that stems from the momentum transfer 
of 2hk per one reflected photon with the wavcnumbcr 
k. It has been been shown that this system may be 
used to generate sub-Poissonian light in the output from 
the cavity O-na]. The moving mirror alters the photon 
statistics by changing the optical path length in a way 
that is proportional to the instantaneous photon number 
inside the cavity. This system may also be used to cre- 
ate highly non-classical states of the cavity field, such as 
Schrodinger cats ||l^,|l3 ; a-nd might even be used to cre- 
ate cat states of the mirror |1^ . In addition, it has been 
shown that such a configuration may be used to perform 
QND measurements of the light field [Q , and to detect 
the decoherence of the mirror, a topic of fundamental 
interest in quantum measurement theory [|l9| . Due to re- 
cent technological developments in optomechanics, this 
area is now becoming experimentally accessible. Dorsel 
et al. have realised optical bistability with this system 
pU, and other experiments, particularly to probe the 



quantum noise, are now in progress |21 2 



In order to create displacements that are large enough 
to be observed, one is tempted to use a mirror having a 
well-defined mechanical resonance with a very high qual- 
ity factor Q. Thus, even when excited with weak white 
noise driven radiation pressure, the mirror can be dis- 
placed by a detectable amount at the mechanical reso- 
nance frequency v. For such a mirror to behave fully 
quantum mechanically one needs to operate at very low 
temperatures since the thermal energy kT very easily ex- 
ceeds hi/. For example, a i'/2tt = 100 kHz resonance is 
already significantly excited at 5 /iK. However, it is not 
necessary to reach the fully quantum domain to observe 
the quantum back action. By simultaneously combining 
a high optical quality factor (ie. by using a high-finesse 
cavity) and a specially designed low mass mirror with 
very high mechanical quality factor one can at typical 
cryogenic temperatures create conditions where the radi- 
ation pressure fluctuations (which are the source of the 
quantum mechanical back-action referred to earlier) ex- 
ceed the effects caused by thermal noise. In this paper we 
discuss considerations for detecting this quantum back- 
action noise. 

There are already a number of publications dealing 
with quantum noise in optical position measurements. 
Our main purpose here is to extend this literature in 
two ways which are important when considering the de- 
tection of the quantum noise. The first is the inclusion 
of the effects of experimental sources of noise, such as 
the classical laser noise and the noise from intracavity 
losses. The second is to perforin a quantum treatment of 
phase-modulation detection, so that the results may be 
compared with those for homodyne detection. While this 
method of phase detection is often used in practice, it has 
not previously been given a quantum mechanical treat- 



ment, which we show is important because previous semi- 
classical treatments have underestimated the shot noise. 
In addition to these main objectives, we also show that 
the standard Brownian motion master equation is not 
adequate to describe the thermal damping of the mirror, 
but that the corrected Brownian motion master equation 
derived by Diosi [E3] rectifies this problem. 

In section || we describe the configuration of the sys- 
tem. In section |l^ we perform a quantum mechanical 
analysis of phase modulation detection. In section IV 
we solve the linearised equations of motion for the cav- 
ity/mirror system, using a non-standard Brownian mo- 
tion master equation which is of the Lindblad form P3| . 
In section M we use this solution to obtain the noise power 
spectral density (which we refer to simply as the spec- 
trum) for a measurement of the phase quadrature using 
phase modulation detection. In the first part of this sec- 
tion we discuss each of the contributions and their respec- 
tive forms. Next we compare the spectrum to that which 
results if the standard (non-Lindblad) Brownian motion 
master equation is used to describe the thermal damping 
of the mirror, and also to that which would have been 
obtained using homodyne detection rather than phase 
modulation detection. Finally we show how the error in 
a measurement of the position of the mirror may be ob- 
tained easily from the spectrum. We evaluate explicitly 
the contribution to this error from various noise sources, 
and plot these as a function of the laser power. Section^ 
concludes. 



II. THE SYSTEM 



The system under consideration consists of a coher- 
ently driven optical cavity with a moving mirror which 
will be treated as a quantum mechanical harmonic oscil- 
lator. The light driving the cavity refiects off the moving 
mirror and therefore fluctuations in the position of the 
mirror register as fluctuations in the light output from 
the cavity. In the limit in which the cavity damping rate 
is much larger than the rate of the dynamics of the mirror 
(characterised by the frequency of oscillation v and the 
thermal damping rate F) the phase fluctuations of the 
output light arc highly correlated with the fluctuations 
of the position of the mirror and constitute a continuous 
position measurement of the mirror Q . 

An experimental realisation will therefore involve a 
continuous phase-quadrature measurement of the light 
output from the cavity to determine the output spectrum 
of the phase-quadrature fluctuations. The nature of the 
detection scheme used to measure the phase quadrature 
is of interest to us, as we shall see that it will effect the 
relationship of the shot noise to the other noise sources 
in the measured signal. Quantum theoretical treatments 
usually assume the use of homodyne detectio n ||7| , p^[ -p^ . 
However this is often not used in practice |p3,p5l. Many 
current experiments use instead phase 
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FIG. 1. The light output from the laser is phase-modulated 
at frequency A in the modulator M, and from there drives the 
cavity. The front mirror of the cavity is fixed, while the back 
mirror is a mechanical harmonic oscillator. The diagonal line 
is a shorthand representation for the arrangement which iso- 
lates the laser from the cavity output, and the light that is 
reflected back off the front mirror. All this light falls upon a 
photo-detector PD, and the photo-detection signal is demod- 
ulated (to pick out the phase quadrature signal) before going 
to a spectrum analyser. 

modulation detection ||2l|| , which was developed by Bjork- 
lund p6| , p7t in 1979. Before we treat the dynamics of the 
cavity field/oscillating mirror system, to determine the 
effect of various noise sources, we will spend some time 
in the next section performing a quantum mechanical 
treatment of phase modulation detection. We will focus 
on this scheme throughout our treatment, and compare 
the results with those for Homodyne detection. A dia- 
gram of the experimental arrangement complete enough 
for the theoretical analysis is given in figure |l|. We note 
that in practice a feedback scheme is used to lock the 
laser to the cavity so as to stabilise the laser frequency. 
For an analysis of the method and an expression for the 
resulting classical phase noise the reader is referred to 
references Q and |2^. We do not need to treat this 
feedback explicitly, however. Its effect may be taken into 
account by setting the value of the classical laser phase 
noise in our analysis to the level it provides. 



III. PHASE MODULATION DETECTION 

The laser which drives the cavity is isolated from the 
cavity output, and the entirety of this output falls upon a 
photo-detector. In order that the photo-detection signal 
contain information regarding the phase quadrature, the 
laser field is modulated at a frequency A, which is chosen 
to be much greater than the natural frequency of the har- 
monic mirror. The sidebands that result from this modu- 
lation are chosen to lie far enough off-resonance with the 
cavity mode that they do not enter the cavity and are 
simply reflected from the front mirror. From there they 
fall upon the photo-detector. The result of this is that 
the output phase quadrature signal appears in the pho- 
todetection signal as a modulation of the amplitude of a 
'carrier' at frequency A. This is then demodulated (by 
multiplying by a sine wave at the modulation frequency 
and time averaging) to pick out the phase quadrature 
signal, and from there the spectrum may be calculated. 



First consider the laser output field, which is essen- 
tially classical; it is a coherent state in which the ampli- 
tude and phase are not completely stable and therefore 
contain some noise. This means that the field from the 
laser actually contains frequencies in a small range about 
its central frequency. The laser field may therefore be 
described by a set of coherent states with frequencies in 
this range. As a result it is possible to perform a uni- 
tary transformation on the mode operators such that the 
amplitude of each of the coherent states is replaced by a 
complex number, and the quantum state of the field is 
simply the vacuum 1 28 3^] . This separates out the classi- 
cal variations in the field from the quantum contribution, 
and allows us to write the output from the laser as 



P + Sain{t) + 6x{t) + iSy{t) 



(1) 



In this expression /3 is the average coherent amplitude of 
the field, which we choose to be real, and /3^ is the pho- 
ton flux. The deviations from this average are given by 
6x{t), being the classical amplitude noise, and Sy{t), be- 
ing the classical phase noise. The quantum noise, which 
may be interpreted as arising from the vacuum quantum 
field, is captured by the correlation function of the field 
operator 6aia{t). Here the subscript refers to the field's 
relation to the cavity, and not the laser. The correlation 
functions of the various noise sources are 

{dai^{t)Sal{t + T)) = S{r) 
{Sal{t)Sai^{t + t)) = 
{Sxit)dx{t + t)) =G^{t) 
{Sy{t)Sy{t + r))=Gy{T), 

where we have left the classical noise sources arbitrary. 
This allows them to be tailored to describe the output 
from any real laser source at a later time. However, we 
will assume that the modulation frequency. A, is choosen 
large enough so that the classical noise is negligible at this 
frequency. This is what is done in practice. The average 
values of the three noise sources, Sai^, Sx and Sy are zero, 
as are all the cross correlations. 

Before the laser field enters the cavity, it passes 
through a phase modulator. This is a classical device 
which modulates the phase of the coherent amplitude of 
the beam, and as such leaves the quantum noise unaf- 
fected. The phase is modulated sinusoidally, the result 
of which is to transform the time dependent coherent 
amphtude, given by /3 + Sx{t) + i6y{t), into P6| , pi| 



{(3 + Sx{t) + tSy{t)) Y. Jn{M)e' 



lAt 



(2) 



n—~oo 



where J„ is the n*'' Bessel function, A is the frequency of 
the modulation, and M is referred to as the modulation 
index, being determined by the amplitude of the sinu- 
soidal modulation. For phase modulation detection, the 
modulation index is typically chosen to be much less than 



unity so that Jq w 1, J±i = zLM/2 = zLe, e <C 1, and 
all other terms vanish. The laser field after modulation 
is then 



{/3 + 6x{t) + tSy{t)){l 



- ee 



-iAt 



ee 



i.At 



) + 5ai„(i). (3) 



Using now the input-output relations of Collet and Gar- 
diner p^ , the field output from the cavity is 



t{t)^-{p + Sx{t)+idy{t))il^ 
-~Sain{t) + y/^{6a{t) + a 



^ -iAt ^ iAt\ 

- ee — ee 



(4) 



in which a{t) = 6a{t) -f a is the operator describing the 
cavity mode, and 7 is the decay constant of the cavity 
due to the input coupling mirror. We are interested in 
the steady state behaviour, and we choose a to be the 
average steady state field strength in the cavity. In ad- 
dition, in order to solve the equations of motion for the 
cavity we will linearise the system about the steady state. 



which requires that {Sa\t)5a{t)) <C 



ar 



The operator 



describing the photo-current from the photo-detector is 



I(t) = flout (i)^ flout (i) 

= &^ + a{6Xout -2Sx) + 2e(3sm{At) ( 6Y, 



2y/jaSy 

/3 



+ {2epysm'{At) 1 + 



2Sx 



(5) 



in which 



and 



5Yont{t)^V^6Y{t)-SYir,{t) 

SX{t) ^ Sa{t) + Sa"! (t) 
SY{t)^-i{6ait)-Sa^t)) 
SX,^{t) ^ 6ai,,{t) + 6alit) 



(6) 



We have also set a = {—f3 + ^/ja), and assumed this to 
be real. To obtain the phase quadrature signal we de- 
modulate, which involves multiplication by a sine wave 
at frequency A, and subsequent averaging over a time, 
T. This time must be long compared to 1/A, but short 
compared to the time scale of the phase quadrature fluc- 
tuations. The signal is therefore given by 



R{t) 



ly sin(AT)/(t + r) dr, 



(7) 



We now must evaluate this to obtain R{t) explicitly in 
terms of the phase quadrature. Writing out the integral, 
and dropping everything which averages to zero (that is, 
which is not passed by the low-pass filtering) we obtain 



R{t) 



SYout{t + T)dT 



+ qi{t) + q2{t), 



where 



6y(t + t) dr 



(8) 



,rW^-('^-^)Re 



T 



92(t)=ff iRe 



ie-^Mt+r)gXinit + T)dT 



~^2A(t+r} ^Yin^t + t) dT 



(9) 
(10) 



In deriving this expression we have assumed that the 
classical laser noise is Note that we choose T to be 
much smaller than the time-scale upon which SY and Sy 
change, so that the integration is essentially equivalent 
to multiplication by T, an effect which is canceled by the 
division by T. However, we should also note that (^Kout 
contains SYin, so that in replacing the first term in R(t) 
by e/3 i5Fout {t) we must remember that this only contains 
the frequency components of SY\n in a bandwidth of 1/T 
around zero frequency. The result of this is that SYout (t) 
is uncorrelated with qi and (72: being the quantum noise 
in the bandwidth 1/T around the frequencies A and 2 A 



respectively. We need to know the correlation functions 
of these noise sources, and whether or not they are cor- 
related with any of the other terms in R{t). It is clear 
that qi and q2 are not correlated over separation times 
greater than 2T. Using Eq.(lO) to evaluate the correla- 
tion function of (72, for example, we have 



{q2{t)q2{t + t)) 



(e/3)' 







for 



<T 



otherwise. 



(11) 



On the time-scale of the fluctuations of SY we can ap- 
proximate this as a delta function, so that gi and 92 (and 
also SYin) are still effectively white noise sources. We may 



therefore write 

R{t) = ef3 SYoutit) + (71 (i) + g2(i) - 2ey^aSyit), (12) 

and the correlation functions of qi and (72 are 

(gi(t)qi(t + r)> = (l/2)(/3 - V^a)2 5(r), 
{q2{t)q2{t + T)) = il/2){ePf6{r). (13) 

The signal therefore contains the phase quadrature of the 
output field, SYoat{t), plus three noise terms. While the 
last term, being the input classical phase noise, is corre- 
lated with SYout{t), qi and 92 are not. Taking the Fourier 
transform of the signal. 



1 r°° 

R{uj) = -j= / R{t)e~"^' dt, 
V 27r J-00 



(14) 



we may write 

R{lj) = 6(3 SYoutH + Vgj(cj) - 2e^a5y{u;). (15) 



2 



=1 



This is the Fourier transform of the signal in the case 
of phase modulation detection. If we were to use ideal 
homodyne detection this would be instead [B6l 



Rh{u;)^Kp{SY,^t{L^)-2Sy{Lu)), 



(16) 



where /3 is the amplitude of the local oscillator and n 
is the reflectivity of the beam splitter used in the ho- 
modyne scheme. Thus, in the case of phase modulation 
detection, there are two white noise sources which do not 
appear in homodyne detection. They stem from the fact 
that the phase quadrature detection method is demod- 
ulating to obtain a signal at a carrier frequency. Be- 
cause the quantum noise is broad band (in particular, it 
is broad compared to the carrier frequency) the demod- 
ulation picks up the quantum noise at A and 2 A. There 
is also a term from the classical phase noise in the side- 
bands. We note that the contribution from the quantum 
noise at 2 A has been omitted from previous semiclassi- 
cal treatments, with the result that the shot noise has 
been underestimated by (1/2) (e/3)^ [g^- For unbalanced 
homodyne detection there will also be an extra contri- 
bution from the noise on the local oscillator, which may 
be suppressed (in the limit of an intense local oscillator) 
with the use of balanced homodyne detection |^^ . 

Returning to Eq.(n3) for the demodulated signal, the 
next step is to solve the equations of motion for the sys- 
tem operators to obtain SY(uj) in terms of the input noise 
sources. We can then readily calculate {R{uj)R{uj')), 
which appears in the form 



{R{cj)R{u;')) ^S{cj)S{u; + cj') 



(17) 



The delta function in to and to' is a result of the stationar- 
ity of R(t), and S{(jj) is the power spectral density, which 
we will refer to from now on simply as the spectrum. This 



is useful because, when divided by 2tt, it gives the aver- 
age of the square of the signal per unit frequency (The 
square of the signal is universally referred to as the power, 
hence the name power spectral density). Since the noise 
has zero mean, the square average is the variance, and 
thus the spectrum provides us with information regard- 
ing the error in the signal due to the noise. The spectrum 
is also a Fourier transform of the autocorrelation func- 
tion [ p3[ . The specific relation, using the definitions we 
have introduced above, is 



Sito) 



(i?(0)i?(r))e- 



^dr. 



(18) 



and as the autocorrelation function has units of s~^, the 
spectrum has units of s^^. To determine the spectrum 
experimentally the phase of the signal is measured for a 
time long compared to the width of the auto-correlation 
function, and the Fourier transform is taken of the result. 
Taking the square modulus of this Fourier transform, and 
dividing by the duration of the measurement obtains a 
good approximation to the theoretical spectrum. We pro- 
ceed now to calculate this spectrum. 



IV. THE DYNAMICS OF THE SYSTEM 

Excluding coupling to reservoirs, the Hamiltonian for 
the combined system of the cavity mode and the mirror 
isfra 



H = hujQa'a 



P 
2m 



1 



-mv'''q — hga* aq 



+n{i [E + ^5x{t) + i^5y{t)] a^ + H.c.} . (19) 

In this equation loo is the frequency of the cavity mode, 
q and p are the position and momentum operators for 
the mirror respectively, m and v are the mass and an- 
gular frequency of the mirror, g — uJa/ L is the coupling 
constant between the cavity mode and the mirror (where 
L is the cavity length), and a is the annihilation opera- 
tor for the mode. The classical driving of the cavity by 
the coherent input field is given by E which has dimen- 
sions of s~^, and is related to the input laser power P 
hy E — y/Pj/Jhujo) = y/jP- The classical laser noise 
appears as noise on this driving term. 

The moving mirror is a macroscopic object at tempera- 
ture T, and as such is subject to thermal noise. While it is 
still common to use the standard Brownian motion mas- 
ter equation (SBMME) |p3| , |40| to model such noise, as it 
works well in many situations, it turns out that it is not 
adequate for our purposes. This is because it generates 
a clearly non-sensical term in the spectrum. As far as we 
know this is the first time that it has been demonstrated 
to fail in the steady state. Discussions regarding the SB- 
MME, and non-Lindblad master equations may found 
in references [p4lB5|. We will return to this point once 
we have calculated the spectrum. We use instead the 



corrected Brownian motion master equation (CBMME) 
derived by Diosi |2^] , to describe the thermal damping of 
the mirror, as this corrects the problems of the SBMME. 
In particular, we use the CBMME in which the cut-off 
frequency of the thermal reservoir is assumed to be much 
smaller than k^T /h. For current experiments k^T /Ti is 
greater than 10 GHz, so this assumption appears reason- 
able, and leads to the simplest Lindblad-form Brownian 
motion master equation. Using this CBMME, and the 
standard master equation for the cavity losses (both in- 
ternal and external) , the quantum Langevin equations of 
motion for the system are given by 



a ^ -j^[a,H] - 



7 + M 



q = --[q,H] + h{r/6mkT)^r]{t) 



a + ^SUinit) + V7*5in(i) (20) 

(21) 
(22) 



p = -j\p,H] - Vp+{2mTkT)-^m, 

in which the correlation functions for the Brownian noise 
sources are 



{v{t)v{t')) = s{t-t'), 



n 



(23) 
(24) 
(25) 



{r,{t)m)^iiV3/2)S{t~t') 



(26) 



In these equations 7 is the decay constant describing 
transmission through the input coupling mirror. All 'in- 
ternal' cavity losses including absorption, scattering and 
loss through the movable mirror are included separately 
via the decay constant /i, and the corresponding vacuum 
fluctuations via the operator bin{t). The effect of me- 
chanical damping and thermal fluctuations of the mirror 
are given by the noise sources £,{t) and ijit) and the me- 
chanical damping constant F. 

We note here that if we were to use the standard Brow- 
nian motion master equation Q.Q, Eqs.(pi|) and ( [2^ ) 
would instead be given by 



■hH] 



(27) 



p^~-[p,H]-Tp+ (2mF/cT)K(0- (28) 

where {C{t)({t')) — S{t — t'). These Langevin equations 
do not preserve the commutation relations of the quan- 
tum mechanical operators, and as a result it is clear that 
the description cannot be entirely correct. 

Calculating the commutators in Eqs.(pO|) to (p2|), we 
obtain 



E 
P 



7 + M 



a + igaq + ^ Sainit) + y/JIbin{t) + ^/jSx{t) +iy/j5y{t) 



+ h{T/6mkT)^r]{t) 

w'^q + hga)a -Tp+ (2mTkT)^£,{t). 



(29) 

(30) 
(31) 



Introducing a cavity detuning Slo (that is, setting the 
cavity resonance frequency in the absence of any cavity 
field to uJc — ujq + Su), and solving these equations for 
the steady state average values we obtain 

^ ' '^ (32, 



{q)ss = 



7 + M 

^9 I 



(33) 



{P)s 







(34) 



where we have set the detuning to Su = g{x)ss to bring 
the cavity on resonance with the driving field in the 
steady state. Linearising the quantum Langevin equa- 
tions about the steady state values, and writing the result 
in terms of the field quadratures, we obtain the following 
linear equations 



/SX\ 

SY 

SQ 
\6PJ 




( V7<5^ir 



(i) + 77Z,5Xb,i„(i) + 

(i) + VM(5rb,in(i) + 
{Tnu/ikT)^r]{t) 
{^TkTl{hv)YH{t) 



2^5x{t) \ 
2^/7<^2/(0 

) 



(35) 



In this set of equations we have scaled the position and 
momentum variables using 



6P = 



mfiv 



{P~ {P)ss)- 



(37) 



2'mi' 



(9 - ('7)ss) 



(36) 



and we have defined x — g{2h/mvy/^, which has units 
of s~^. The quadratures for the input noise due to intra- 
cavity losses are given by 



SYb^in = -i{bin - b}^ 



(38) 
(39) 



Without loss of generality we have chosen the input field 
amplitude to be real (Ini[/3] = 0), so that the input phase 
quadrature is given by Yin. We now solve the dynam- 
ics ( p5[ ) in the frequency domain in order to obtain the 
spectrum directly from the solution. To switch to the 
frequency domain we Fourier transform all operators and 
noise sources. In particular we have, for example 



da{Lu) 
SaHiu) 



1 



/27r 
1 



Sa{t) 



j*'^* dU 



/2n J- 



6a\t)e''^' dt 



[5a(-c.)]t. 



(40) 
(41) 



Rearranging the transformed equations, the solution is 
given by 

iSX{uj), SY{uj), 6Q{uj), 5P{uj)f = M{uj)n{uj) (42) 

where n(a;) is the vector of transformed noise sources. 
If we write the matrix elements of M{lo) as Mij{cu) — 
mij(uj)/ D(uj), then 

D{u;) = ((7 + ^)/2 - iujf {v^ -u? - iVu}) (43) 

and the non-zero mij are given by 



{{^^ + ^)l2-iu:){i?-u? 



mil 



"fn^i = X^a^t^ 



iTll}) 



m22, 



"123 = xa(r - iuj) ((7 + m)/2 - iuj) , 
m24 = Xo:v ((7 + Ai)/2 - i^) = "^31, 
W33 = (r - iuj) ((7 + /i)/2 - iujf , 

W34 = i^ ((7 + m)/2 - i(^f = -TO43, 
"'-41 = -ixauj ((7 + /i)/2 - iuj) , 
TO44 = — io; ((7 + /i)/2 — itj) . 



(44) 



We have now solved the equations of motion for the sys- 
tem in frequency space. The spectra of the system vari- 
ables may now be calculated in terms of the input noise 



sources. Using the input-output relations, which give the 
output field in terms of the system variables and the in- 
put noise sources, the spectra of the output field, and 
hence of the measured signal, may be obtained. Note 
that quantum mechanics plays no role in the solution of 
the motion of the system. The linear equations of motion 
may as well be equations for classical variables. The only 
part that quantum mechanics plays in determining the 
spectra of the system variables is that some of the input 
noise sources are quantum mechanical. That is, their cor- 
relation functions are determined by quantum mechanics. 
In fact, if all the noise sources had purely classical cor- 
relation functions, then the SBMME Langevin equations 
would not lead to any problems, as they are perfectly 
correct as equations of motion for a classical system. 



V. THE POWER SPECTRAL DENSITY 



To calculate the spectrum of the signal, we require the 
correlation functions of the input noise sources. To reit- 
erate, these are 

{6Xin{u;)6Xin{Lu')) - {6Yi^iLu)6YiniLo')) = 6(0; + u'), 
{SXU^)SYi^{Lu')) - ^{SYiniij)SXi^iLu')) = iSiu + u'), 

and similarly for (5Xb,in(w) and JYb,in(i^)- The correla- 
tion functions for the classical laser noise, and thermal 
noise sources are 

{Sx{uj)Sx{uj')) = Gx{uj)S{u; + uj'), 
{Sy{uj)6y{uj'))^Gy{uj)S{LU + uj'), 

ivimt')) = -imvit')) = ^{V3/2)Sit t'). (45) 

After some calculation we obtain the spectrum of the 
signal for phase modulation detection as 



{epy 
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+ ^Gy[LO) 
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+ 7(7 + Ai + 47G^(w)) 



UJ 



(7 + A^)' \(^l+Ef+uj\ 



-^{xaYT[Av'Ts + -{T'+u:')T^ 






;4H)^+^2 



\D{uj) 



(46) 



where 



\D{uj)\^ = [(7 + m) V4 + -^] ' [(-^ - uj^f + rV] , 



(47) 



and Ts is a dimensionless scaled temperature given by 
2s = (k^ / {hi>))T . This phase-fluctuation spectrum may 
be thought of as arising in the following way. The me- 
chanical harmonic oscillator, which is the moving mirror. 



is driven by various noise sources, both quantum me- 
chanical and classical in origin, and the resulting position 
fluctuations of the mirror are seen as fluctuations in the 
phase of the light output from the cavity. 

Let us examine the origin of the various terms in the 
spectrum in turn. The first two terms, which appear in 
the first set of square brackets, are independent of the 
frequency, and are the contribution from the (quantum 
mechanical) shot noise of the light. The first term has the 
factor of three (rather than a factor of two which would 
be the case for homodyne detection) due to the contri- 
bution from q2{t). The second term is the contribution 
from qi (t) . 

The next three terms, which multiply the second set of 
square brackets, are the back-action of the light on the 
position of the mirror, noise from internal cavity losses, 
and the classical amplitude noise on the laser, respec- 
tively. Note that the only distinction between the back- 
action and the internal losses is that former is propor- 
tional to the loss rate due to the front mirror, and the 
latter is proportional to the internal loss rate. It is eas- 
ily seen that these noise sources should have the same 
effect upon the position of the mirror: the back-action 
is due to the random way in which photons bounce off 
the mirror, whereas the internal losses are due to the 
similarly random way in which photons are absorbed by 
the mirror, (or anything else in the cavity). The am- 
plitude fluctuations of the laser also affect the mirror in 
the same manner, but since these fluctuations are not 
white noise (as is the case with the quantum noise which 
comes from the photon 'collisions'), the response func- 
tion of the mirror is multiplied by the spectrum of the 
amplitude fluctuations. 

The term which appears in the third set of square 
brackets is due to the classical phase fluctuations of the 
laser. Clearly this has quite a different form from that 
due to the quantum noise and the classical amplitude 
fluctuations. In particular, it is not dependent upon the 
coupling constant, g, because it is derived more or less 
directly from the input phase noise. Conversely, the noise 
that derives from the amplitude fluctuations has its ori- 
gin from the fact that the amplitude fluctuations first 
drive the mirror, and it is the resulting position fiuctu- 
ations which cause the phase fiuctuations in the output. 
The classical phase noise term includes a contribution 
from the laser phase noise refiected from the cavity (that 
is, the term given explicitly in Eq.(p^), and a contribu- 
tion from the phase noise on the light which has passed 
through the cavity (being a part of l^ut)- 

The final two terms, which multiply the fourth set of 
square brackets, are due to the thermal fiuctuations of 
the mirror. Note that these terms are only valid in the 
region in which fee 7" ^ hv. 

Finally we note that we do not see squeezing in the 
spectrum of phase quadrature fluctuations. This is be- 
cause squeezing is produced when the cavity detuning is 
chosen so that the steady state detuning is non-zero Q . 
We have chosen to set the steady state detuning to zero 



in this treatment as we are not concerned here with re- 
ducing the quantum noise. 

In what follows we examine various aspects of the spec- 
trum which are of particular interest. Before discussing 
considerations for detecting the back-action noise, we 
compare the spectrum with that which would have been 
obtained using the SBMME, and for that which would 
result from the use of homodyne detection. We then 
write the spectrum at resonance as a function of the laser 
power, and plot this for current experimental parame- 
ters. So far we have been considering the noise power 
spectrum, and have made no particular reference to the 
limit this implies fo r a m easurement of the position of the 
mirror. In section |VC] we show how the spectrum tells 
us the limit to the accuracy of position measurement in 
the presence of the noise sources. 



A. Comparison with the Standard Treatment of 
Brownian Motion 



To obtain the spectrum we have used the corrected 
Brownian motion master equation p3| . This is essential 
because the spectrum which results from the standard 
Brownian motion master equation contains a term which 
is assymetric in uj, and therefore clearly incorrect. In 
particular, to obtain the spectrum given by the SBMME 
from that given by the CBMME, the term proportional 
to Tg must be replaced by 



2u!jTx a V 



((^)'+-^) 



\D{lo) 



(48) 



That the spectrum must be symmetric in oj follows read- 
ily from the stationarity of the output field, and the fact 
that the output field commutes with itself at different 
times. In particular, the stationarity of the output field 
means that the correlation function of the signal only 
depends on the time difference, so that 



{R{t)R{t + T)) =G{t). 



(49) 



As the output field commutes with itself at different 
times, R commutes with itself at different times, and we 
have 

G{-t) = {R{t)R{t - t)) = {R{t - T)R{t)) = G{t). (50) 

The correlation function is therefore symmetric in t. As 
the spectrum is the Fourier transform of the correlation 
function, it follows from the properties of the Fourier 
transform that the spectrum is symmetric in uj. 

It was shown in Ref. ||3g] that for realistic systems at 
high temperatures the SBMME has a stationary density 
matrix which is positive. The non-Lindblad nature of 
the master equation appears only to cause problems at 
short times. In our problem we are calculating spectra 



at steady-state so it might seem surprising that the non- 
Lindbad nature does cause problems for us. On reflec- 
tion, however, this is not surprising. The spectra we cal- 
culate are for continuously measured quantities. Making 
such measurements continuously reprepares the system 
in a conditioned state which is different from the station- 
ary state. Thus if one is observing the system then it is 
never really at steady state and the "initial slip" problem 
of Ref. iQ never goes away. 

Diosi's corrected Brownian motion master equation re- 
moves the term asymmetric in lu by adding a noise source 
to the position (see Eq.(|2^)) which is correlated with the 
noise source for the momentum. In doing so it produces 
an additional term in the spectrum proportional to l/T, 
an effect which, it should be noted, is independent of the 
phase detection scheme. For temperatures (and frequen- 
cies) for which this new term is much smaller than the 
standard term, which is proportional to T, this new term 
can be neglected. However, the question of observing this 
term experimentally is a very interesting one, because it 
would allow the CBMME to be tested. Comparing the 
new term with the term proportional to T wc find that 
the new term begins to dominate when 



T< 



12kp 



Vf^ 



(51) 



For temperatures of the order of a few Kelvin, the addi- 
tional term therefore becomes apparent in the spectrum 
at frequencies of a few Gigahertz. Note that for such high 
frequencies phase modulation may no longer be practical 
however, owing to the fact that A must be much larger 
than the frequency range of the signal. In that case the 
use of alternative phase detection schemes would be re- 
quired 



B. Comparison with Homodyne Detection 

Let us now briefly compare the spectrum derived above 
for phase modulation detection to that which would be 
obtained with homodyne detection. Firstly, if homodyne 
detection had been used, the overall scaling of the spec- 
trum would be different, as it would be proportional to 
the strength of the local oscillator. Thus the factor of 
l/(e/3)^ would be replaced by 1/(/3k)^, in which /3 and 
K are as defined in Eq. (|l(^). This overall factor aside, 
two terms in the spectrum would change. The shot noise 
component would be reduced to unity, and the classical 
phase noise contribution would become 



AGyico) 






(52) 



C. The Error in a Measurement of Position 

So far we have been considering the noise spectrum of 
the phase quadrature, as this is what is actually mea- 
sured. In this section we show how the error in a mea- 
surement of the position of the mirror may be obtained 
in a simple manner from the spectrum, Eq.(|46|), and give 
an example by calculating it explicitly for some of the 
terms. As explained above, the reason for performing 
the phase measurement is that it constitutes essentially 
a measurement of the position of the mirror. 

We can choose to measure the amplitude of position 
oscillations at any frequency, but for the purposes of dis- 
cussion, a measurement of a constant displacement is the 
simplest. First we must see how the position of the mirror 
appears in the signal, which is the phase quadrature mea- 
surement (that is, convert from the units of the signal into 
units of the position fluctuations) . This is easily done by 
calculating the contribution to the spectrum of the posi- 
tion fluctuations due to one of the noise sources (for the 
sake of definiteness we will take the thermal noise) , and 
comparing this to the equivalent term in the spectrum 
of the signal. This gives us the correct scaling. Perform- 
ing this calculation, we find that the spectrum of position 
fluctuations of the mirror due to thermal noise is given by 
the thermal term in the spectrum (Eq.(pq)), multiplied 
by the factor 



7 + M 



2m[e0)'^^i'x'^oi 



2rv2 



(53) 



From this we see that the scaling factor is frequency de- 
pendent. This means, that the spectrum of the position 
fluctuations is somewhat different from the spectrum of 
the resulting phase quadrature fluctuations. For the mea- 
surement of the phase to correspond to a true measure- 
ment of the position the two spectra should be the same. 
This is true to a good approximation when 7 is much 
larger than the range of uj over which the spectrum of po- 
sition fluctuations is non-zero, and this is why the scheme 
can be said to constitute a measurement of position when 

7> i^,r. 

In performing a measurement of a constant displace- 
ment of the mirror (achieved by some constant external 
force), the signal (after scaling appropriately so that it 
corresponds to position rather than photocurrent) is inte- 
grated over a time t^. The best estimate of the displace- 
ment is this integrated signal divided by the measurement 
time. The error. Ax, in the case that the measurement 
time is much greater than the correlation time of the 
noise, is given by 



Aa:2(0)= / (i?,(0)i?,(r))dr/r,n = 5,(0)/7 



(54) 



In this equation R^ and Sx are the appropriately scaled 
signal and spectrum. To calculate the error in the mea- 
surement of a constant displacement, all we have to do. 



therefore, is to scale the spectrum using the expression 
Eg. (p3[) , evaluate this at zero frequency, and divide by 
the measurement time. In general, the spectrum eval- 
uated at a given frequency, once divided by the mea- 
surement time, gives the error in a measurement of the 
amplitude of oscillations at that frequency. We calcu- 
late now the contribution to the error in a measurement 
at zero frequency and at the mirror resonance frequency, 
from the shot noise, thermal, and quantum back-action 
noise. In the following we write the expressions in terms 
of the parameters usually used by experimentalists: the 
laser power, P, cavity finesse JF, and the quality factor 
for the mirror oscillator, Q — v/T. We chose the cavity 
to be impedance matched, since this is usually the case 
in practice. This means that the decay rate due to the 
input coupler, 7, is chosen equal to the internal cavity 
decay rate, /z. The total decay rate of the cavity is there- 
fore 27, so that the finesse is given by J-" = 'nc/{2L^). 
We also assume that 7^1/, which is certainly true in 
current experiments. Performing the calculation we find 
that the contribution due to the shot noise is the same 
at all frequencies, and is given by 



Error in a position measurement vs. Laser Power 
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(55) 



The contribution from the quantum back action for a 
measurement of a constant displacement is 



Ax?,^(0) 



2 T ,4 



BAV'^y — 9 I 9 / 



(56) 



and for a measurement at the resonance frequency v it is 
A^ba(^) = Q^A^baC^)- Note that since ^i = j, the con- 
tribution from the internal cavity losses is also given by 
this expression. In a sense, the internal cavity loss noise 
can also be regarded as a back-action term, although the 
back action is from a measurement process due to the in- 
teraction with an environment that is not being observed. 
The total error which can be said to arise from the ran- 
dom 'photon impacts' on the mirror (in the absence of 
classical laser noise) is the sum of the back-action and 
internal loss noise, and is therefore given by 



a4n(o) 



\c' J 



1 



T^P 



The contribution from the thermal noise is 



A4h(0) = 



2kBT 



QmvkBTQ^Tn 



for a constant displacement, and is 

'2kBQT\ f h^c 



Ax^i^ii^) 



\6miykBTT„ 



(57) 



(58) 



(59) 



for an oscillation at the mirror frequency. In obtaining 
the second term in this last expression we have also used 




FIG. 2. The error in a measurement of the position of the 
mirror, for a measurement time of r^ ~ 300s. The dashed 
curve corresponds to a measurement of a constant displace- 
ment, and the dot-dash curve to a measurement at the mirror 
resonance frequency. The combined contribution of the quan- 
tum back-action and internal cavity loss noise to both curves 
is the sloping section to the right. The quantum back-action 
noise may be observed at reasonable laser power levels. For 
parameters see text. 

i/ 2> r. The contribution from the other noise sources 
may also be readily evaluated from the terms in the spec- 
trum Eq.(|6|). 

Let us examine the total error in a position measure- 
ment resulting from these four contributions (shot-noise, 
back-action, internal losses, and thermal noise) for state- 
of-the-art experimental parameters. Reasonable values 
for such parameters are as follows pi[. The laser fre- 
quency is ti>o = 27r X 2.82 x 10^^ rad s^ (assuming a 
Nd:YAG laser with a wavelength of 1064 nm), the cavity 
length is L = 1 cm, the mass of the oscillating mirror is 
m = 10~^ kg, and the resonant frequency of the mirror is 
v = 2nx2x 10** rad s~ . The quality factor of the mirror 
is four million, which gives F « 3 x 10^^ s^^. With these 
parameters for the cavity we have x = 2.29 x 10~* s~^. 
The cavity damping rate through the front mirror is 
7 = 4.7 X 10^ s^^, and we assume impedance matching so 
that /i = 7. The cavity may be cooled to a temperature 
of T = 4.2 K, so that Tg = kBT/{hv) = 4.37 x 10^ which 
is certainly in the high temperature regime (Tg ^ 1). 
The Diosi term (of order Tg"^ at resonance) is thus to- 
taly negligible. 

In figure g we plot the position measurement error as 
a function of the laser power, both for the measurement 
of a constant displacement, and for a displacement at 
the mirror resonance frequency. The expressions for the 
measurement error derived above are valid in the limit 
where the measurement time is much greater than the 
correlation time of the noise. As the cavity-mirror sys- 
tem is driven by white noise, this correlation time is given 
approximately by the longest decay time of the system. 
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In our case this is the decay time of the moving mirror, 
given by l/F « 30 s. In view of this we have chosen a 
measurement time of 300 seconds (5 minutes) for the plot 
in figure y. 

The uncertainty due to the shot noise falls off with laser 
power, while that due to thermal noise is independent of 
laser power, and that due to the quantum back-action in- 
creases with laser power. These results are already well 
known. The thermal and back-action contributions are 
much greater at the resonance frequency of the mirror, 
due to the high mechanical Q factor. The optimal regime 
for detecting the quantum back-action noise is at reso- 
nance, as the absolute magnitude of this noise is largest 
in this case. Reasonable experimental values for laser 
power lie between the solid lines, where the increase in 
noise due to the back-action is visible. However, our 
analysis of the spectrum shows us that the full situa- 
tion is more complicated. We have shown that the noise 
due to internal cavity losses and the classical laser am- 
plitude noise have the same dependence on frequency as 
the quantum back-action. In order to reach the back- 
action dominated regime, the laser amplitude noise must 
be at the shot noise level, and the frequency noise must 
be extremely low. 



dent on the phase measurement scheme. In particular, 
the result for phase modulation detection, which is com- 
monly used in experiments of this kind, is not the same as 
that for homodyne detection. We have found that while 
the signature of the classical phase noise is quite differ- 
ent for that of the quantum-back action, the noise due 
to intracavity losses and classical amplitude noise has a 
very similar signature to the back-action. As far as the 
parameters of the cavity and oscillating mirror are con- 
cerned, realisable experiments are beginning to fall in the 
region where the quantum back-action may be observed. 
In our treatment of the system we have shown that 
the standard quantum Brownian motion master equation 
produces a clearly spurious term in the steady state noise 
spectrum for the phase quadrature measurement. We 
have shown that the corrected Brownian motion master 
equation, derived by Diosi, corrects this error. However, 
it also produces a new term in the spectrum which is 
small for present experimental systems. Testing for the 
existence of this term poses an experimental challenge 
that might be met using miniature, high frequency oscil- 
lators and ultra-low temperatures. 
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